Update of explore-JULES-ES-1p0.Rmd that uses the new packages and functions
We build and test gaussian process emulators, perform a sensitivity analyis, and constrain the input space of Jules.
Andy thinks there might be mileage in analysing the atmospheric growth, which is here. /home/h01/hadaw/Research/ES_PPE/Oct20
At the moment, this vignette is hampered by the fact that emulators are failing on a few of the outputs which represent change over the historical period. The emulator is fine for predicting absolute values in the modern period.
Andy would like to see timeseries of: cVeg, cSoil and nbp, npp in GtC and GtC/yr.
Load libraries, functions and data.
# Load helper functions
knitr::opts_chunk$set(fig.path = "figs/", echo = FALSE, message = FALSE, warnings = FALSE)
# load some helper functions
source('~/brazilCSSP/code/brazil_cssp/per_pft.R') # eventually, move the relevant functions
source('explore-JULES-ES-1p0_PPE_functions.R')
Section 2.5 in Friedlingstein et al. describes how the land carbon sink is estimated.
For each ensemble member and each variable, calculate the change from the 20 years at the start of the run, to the twenty years at the end of the run.
Initial clean of data set, removing variables that don’t change, and removing NAs (models that didn’t run).
[1] TRUE
(“probability of run failure” Could be an interesting project - logit transformation for probability of run failure? Some other ML like SVM?
There are clear run failure thresholds in the parameters rootd_ft_io and lai_max, and quite strong visual indications that a_wl_io and bio_hum_cn matter.
NEXT, histograms of failure with the LHS value
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0069/stats/JULES-ES-1p0_P0069_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0086/stats/JULES-ES-1p0_P0086_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0144/stats/JULES-ES-1p0_P0144_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0165/stats/JULES-ES-1p0_P0165_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0168/stats/JULES-ES-1p0_P0168_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0183/stats/JULES-ES-1p0_P0183_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0185/stats/JULES-ES-1p0_P0185_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0214/stats/JULES-ES-1p0_P0214_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0216/stats/JULES-ES-1p0_P0216_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0235/stats/JULES-ES-1p0_P0235_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0238/stats/JULES-ES-1p0_P0238_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0251/stats/JULES-ES-1p0_P0251_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0259/stats/JULES-ES-1p0_P0259_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0285/stats/JULES-ES-1p0_P0285_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0357/stats/JULES-ES-1p0_P0357_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0376/stats/JULES-ES-1p0_P0376_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0387/stats/JULES-ES-1p0_P0387_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0394/stats/JULES-ES-1p0_P0394_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0405/stats/JULES-ES-1p0_P0405_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0424/stats/JULES-ES-1p0_P0424_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0455/stats/JULES-ES-1p0_P0455_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
Error in R_nc4_open: No such file or directory
Error in nc_open(paste0(fn)) :
Error in nc_open trying to open file /project/carbon_ppe/JULES-ES-1p0_PPE/P0457/stats/JULES-ES-1p0_P0457_Annual_global.nc
Error in R_nc4_inq_varndims: NetCDF: Not a valid ID
Error in ncvar_ndims(ncid, varid) : error returned from C call
Error in R_nc4_close: NetCDF: Not a valid ID
First, use mean NPP as an example. How does NPP respond to each parameter? NAs are removed, but zero values are still included.
Some outputs (e.g fLuc, fHarvest) have an almost perfect 1:1 relationship between modern value and change, some (nbp, npp, treeFrac) quite or moderately strong, and some (csoil, cveg) very weak or non-existant.
It appears that this ensemble is less “clear cut” in having an output that clearly distinguishes between “failed” (or close to it), and “not failed”.
Having said that, having an F0 over a threshold seems to kill the carbon cycle, as before. Here, we’ve set a threshold of 0.9 (on the normalised scale) for F0, and we remove members of the ensemble with a larger F0 than that when we build emulators.
The level 1 constraint removes any input with F0 greater than 0.9 (normalised), which removes many of the zero-carbon-cycle members up front. There are 424 ensmble members remaining.
This does constraint sequentially, which may not be a good idea.
[1] 2.410932e-09
# Plot the regular km emulator. Doesn’t look great.
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 2 2 2 2 1.998154 1.997409 2 2 2 2 2 2 2 2 2 1.999201 2 1.996693 2 2 2 2 2 2 1.995643 1.998832 2 2 2 2
- best initial criterion value(s) : -1842.508
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1842.5 |proj g|= 1.8652
At iterate 1 f = 1841.3 |proj g|= 1.4365
At iterate 2 f = 1839.3 |proj g|= 1.8697
At iterate 3 f = 1839 |proj g|= 1.8947
At iterate 4 f = 1839 |proj g|= 1.2143
At iterate 5 f = 1838.6 |proj g|= 1.1891
At iterate 6 f = 1838.3 |proj g|= 1.9015
At iterate 7 f = 1837.6 |proj g|= 1.871
At iterate 8 f = 1837 |proj g|= 1.8711
At iterate 9 f = 1834.9 |proj g|= 1.7796
At iterate 10 f = 1833 |proj g|= 1.8503
At iterate 11 f = 1831.1 |proj g|= 1.8644
At iterate 12 f = 1830.9 |proj g|= 1.8531
At iterate 13 f = 1829.9 |proj g|= 1.8044
At iterate 14 f = 1828.8 |proj g|= 1.8322
At iterate 15 f = 1827.7 |proj g|= 1.8111
At iterate 16 f = 1825.8 |proj g|= 1.8566
At iterate 17 f = 1822.1 |proj g|= 1.8783
At iterate 18 f = 1817.9 |proj g|= 1.8192
At iterate 19 f = 1810 |proj g|= 1.8972
At iterate 20 f = 1807.7 |proj g|= 1.3523
At iterate 21 f = 1806.7 |proj g|= 1.3328
At iterate 22 f = 1804.9 |proj g|= 1.8829
At iterate 23 f = 1803.2 |proj g|= 1.6989
At iterate 24 f = 1802.4 |proj g|= 1.7693
At iterate 25 f = 1801.2 |proj g|= 1.8451
At iterate 26 f = 1799.8 |proj g|= 1.8602
At iterate 27 f = 1799.5 |proj g|= 1.6174
At iterate 28 f = 1799.2 |proj g|= 1.5947
At iterate 29 f = 1798.8 |proj g|= 1.5503
At iterate 30 f = 1796.6 |proj g|= 1.8546
At iterate 31 f = 1795.7 |proj g|= 1.586
At iterate 32 f = 1795.4 |proj g|= 1.583
At iterate 33 f = 1795.3 |proj g|= 1.6047
At iterate 34 f = 1793.7 |proj g|= 1.5063
At iterate 35 f = 1793.4 |proj g|= 1.4599
At iterate 36 f = 1793.2 |proj g|= 1.8832
At iterate 37 f = 1793 |proj g|= 1.87
At iterate 38 f = 1792.9 |proj g|= 1.4333
At iterate 39 f = 1792.4 |proj g|= 1.6437
At iterate 40 f = 1792.1 |proj g|= 1.3598
At iterate 41 f = 1791.8 |proj g|= 1.8392
At iterate 42 f = 1791.2 |proj g|= 1.8717
At iterate 43 f = 1790.7 |proj g|= 1.8739
At iterate 44 f = 1790.5 |proj g|= 1.0022
At iterate 45 f = 1790.4 |proj g|= 1.1185
At iterate 46 f = 1790.3 |proj g|= 1.8759
At iterate 47 f = 1790.3 |proj g|= 1.875
At iterate 48 f = 1790.2 |proj g|= 1.8727
At iterate 49 f = 1790.1 |proj g|= 1.802
At iterate 50 f = 1790.1 |proj g|= 0.81547
At iterate 51 f = 1790.1 |proj g|= 0.8373
At iterate 52 f = 1790 |proj g|= 1.8728
At iterate 53 f = 1789.8 |proj g|= 0.93276
At iterate 54 f = 1789.8 |proj g|= 1.8773
At iterate 55 f = 1789.7 |proj g|= 1.8754
At iterate 56 f = 1789.7 |proj g|= 1.5203
At iterate 57 f = 1789.6 |proj g|= 1.2658
At iterate 58 f = 1789.6 |proj g|= 0.92604
At iterate 59 f = 1789.2 |proj g|= 1.876
At iterate 60 f = 1789.2 |proj g|= 1.873
At iterate 61 f = 1789.1 |proj g|= 0.80926
At iterate 62 f = 1789.1 |proj g|= 0.63617
At iterate 63 f = 1789.1 |proj g|= 1.7888
At iterate 64 f = 1789 |proj g|= 1.873
At iterate 65 f = 1789 |proj g|= 1.4746
At iterate 66 f = 1789 |proj g|= 0.71729
At iterate 67 f = 1789 |proj g|= 0.71538
At iterate 68 f = 1789 |proj g|= 0.70747
At iterate 69 f = 1789 |proj g|= 0.63382
At iterate 70 f = 1788.9 |proj g|= 0.74475
At iterate 71 f = 1788.9 |proj g|= 1.7914
At iterate 72 f = 1788.9 |proj g|= 0.95343
At iterate 73 f = 1788.8 |proj g|= 0.85278
At iterate 74 f = 1788.8 |proj g|= 0.74093
At iterate 75 f = 1788.8 |proj g|= 0.90557
At iterate 76 f = 1788.8 |proj g|= 0.29949
At iterate 77 f = 1788.8 |proj g|= 0.28806
At iterate 78 f = 1788.8 |proj g|= 0.27628
At iterate 79 f = 1788.8 |proj g|= 0.55663
At iterate 80 f = 1788.8 |proj g|= 1.7132
At iterate 81 f = 1788.8 |proj g|= 1.8713
At iterate 82 f = 1788.8 |proj g|= 1.8716
At iterate 83 f = 1788.8 |proj g|= 1.8729
At iterate 84 f = 1788.8 |proj g|= 0.83666
At iterate 85 f = 1788.8 |proj g|= 0.30137
At iterate 86 f = 1788.7 |proj g|= 0.31882
At iterate 87 f = 1788.7 |proj g|= 0.2255
At iterate 88 f = 1788.7 |proj g|= 0.42148
At iterate 89 f = 1788.7 |proj g|= 0.59752
At iterate 90 f = 1788.7 |proj g|= 0.51451
At iterate 91 f = 1788.7 |proj g|= 1.4762
At iterate 92 f = 1788.7 |proj g|= 0.69297
At iterate 93 f = 1788.7 |proj g|= 0.53596
At iterate 94 f = 1788.7 |proj g|= 0.5359
At iterate 95 f = 1788.7 |proj g|= 0.55123
At iterate 96 f = 1788.7 |proj g|= 0.37871
At iterate 97 f = 1788.7 |proj g|= 0.31548
At iterate 98 f = 1788.7 |proj g|= 0.23863
At iterate 99 f = 1788.7 |proj g|= 0.13991
At iterate 100 f = 1788.7 |proj g|= 0.13728
At iterate 101 f = 1788.7 |proj g|= 0.14306
final value 1788.683595
stopped after 101 iterations
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 2 1.995574 2 1.998154 1.790352 2 2 2 1.998084 2 2 2 2 2 1.999201 2 1.996693 1.998361 2 2 1.995738 2 1.993455 1.995643 1.998832 2 2 2 2
- best initial criterion value(s) : -6192.077
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 6192.1 |proj g|= 1.7779
At iterate 1 f = 6097.2 |proj g|= 1.9547
At iterate 2 f = 6092.5 |proj g|= 1.9408
At iterate 3 f = 6091.7 |proj g|= 1.94
At iterate 4 f = 6089.5 |proj g|= 1.9104
At iterate 5 f = 6081.3 |proj g|= 1.9157
At iterate 6 f = 6077.6 |proj g|= 1.9116
At iterate 7 f = 6074.3 |proj g|= 1.907
At iterate 8 f = 6073.6 |proj g|= 1.9049
At iterate 9 f = 6071.5 |proj g|= 1.8255
At iterate 10 f = 6066 |proj g|= 1.7837
At iterate 11 f = 6065.2 |proj g|= 1.7673
At iterate 12 f = 6063.4 |proj g|= 1.752
At iterate 13 f = 6062.6 |proj g|= 1.9208
At iterate 14 f = 6062.1 |proj g|= 1.7448
At iterate 15 f = 6062 |proj g|= 1.7437
At iterate 16 f = 6061.8 |proj g|= 1.7402
At iterate 17 f = 6061 |proj g|= 1.7269
At iterate 18 f = 6060.1 |proj g|= 1.7134
ys=-1.254e-02 -gs= 8.259e-01, BFGS update SKIPPED
At iterate 19 f = 6059.6 |proj g|= 1.7043
ys=-2.746e-03 -gs= 5.566e-01, BFGS update SKIPPED
At iterate 20 f = 6059.3 |proj g|= 1.7005
At iterate 21 f = 6056.3 |proj g|= 1.9169
At iterate 22 f = 6055.7 |proj g|= 1.9087
At iterate 23 f = 6055.4 |proj g|= 1.8996
At iterate 24 f = 6055.1 |proj g|= 1.6244
At iterate 25 f = 6054.9 |proj g|= 1.8958
At iterate 26 f = 6054.6 |proj g|= 1.6078
At iterate 27 f = 6054.3 |proj g|= 1.5928
At iterate 28 f = 6053.9 |proj g|= 1.8949
At iterate 29 f = 6053.8 |proj g|= 1.8942
At iterate 30 f = 6053.7 |proj g|= 1.5209
At iterate 31 f = 6053.7 |proj g|= 1.516
At iterate 32 f = 6053.3 |proj g|= 1.8931
At iterate 33 f = 6052.9 |proj g|= 1.3806
At iterate 34 f = 6052.8 |proj g|= 1.3433
At iterate 35 f = 6052 |proj g|= 1.8927
At iterate 36 f = 6051.9 |proj g|= 1.3527
At iterate 37 f = 6051.9 |proj g|= 1.0531
At iterate 38 f = 6051.8 |proj g|= 0.96661
At iterate 39 f = 6051.8 |proj g|= 0.97303
At iterate 40 f = 6051.8 |proj g|= 0.62567
At iterate 41 f = 6051.7 |proj g|= 0.68667
At iterate 42 f = 6051.7 |proj g|= 1.0212
At iterate 43 f = 6051.7 |proj g|= 0.36663
At iterate 44 f = 6051.7 |proj g|= 0.20017
At iterate 45 f = 6051.7 |proj g|= 0.20136
At iterate 46 f = 6051.7 |proj g|= 0.1629
At iterate 47 f = 6051.7 |proj g|= 0.21978
At iterate 48 f = 6051.7 |proj g|= 0.23363
At iterate 49 f = 6051.7 |proj g|= 0.17448
At iterate 50 f = 6051.7 |proj g|= 0.091728
At iterate 51 f = 6051.7 |proj g|= 0.097416
At iterate 52 f = 6051.7 |proj g|= 0.12446
At iterate 53 f = 6051.7 |proj g|= 0.07705
At iterate 54 f = 6051.7 |proj g|= 0.11677
At iterate 55 f = 6051.7 |proj g|= 0.20318
At iterate 56 f = 6051.7 |proj g|= 0.094594
At iterate 57 f = 6051.7 |proj g|= 0.11668
At iterate 58 f = 6051.7 |proj g|= 0.03121
At iterate 59 f = 6051.7 |proj g|= 0.055344
At iterate 60 f = 6051.7 |proj g|= 0.11291
At iterate 61 f = 6051.7 |proj g|= 0.066099
At iterate 62 f = 6051.7 |proj g|= 0.078269
iterations 62
function evaluations 76
segments explored during Cauchy searches 76
BFGS updates skipped 2
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0782694
final function value 6051.72
F = 6051.72
final value 6051.720344
converged
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 2 2 2 2 1.998154 1.997409 2 2 2 2 2 2 2 2 2 1.999201 2 1.996693 2 2 2 2 2 2 1.995643 1.998832 2 2 2 2
- best initial criterion value(s) : -1844.124 -1844.169 -1844.496 -1844.519
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1844.1 |proj g|= 1.9783
At iterate 1 f = 1827.3 |proj g|= 1.5081
At iterate 2 f = 1810.4 |proj g|= 1.9257
At iterate 3 f = 1810.1 |proj g|= 1.2762
At iterate 4 f = 1808.6 |proj g|= 1.9155
At iterate 5 f = 1803.8 |proj g|= 1.9384
At iterate 6 f = 1802.5 |proj g|= 1.9352
At iterate 7 f = 1801.9 |proj g|= 1.9195
At iterate 8 f = 1801.3 |proj g|= 1.5673
At iterate 9 f = 1801 |proj g|= 1.5746
At iterate 10 f = 1800.3 |proj g|= 1.5734
At iterate 11 f = 1800 |proj g|= 1.5698
At iterate 12 f = 1799.2 |proj g|= 1.9053
At iterate 13 f = 1799 |proj g|= 1.9022
At iterate 14 f = 1798.9 |proj g|= 1.2078
At iterate 15 f = 1798.9 |proj g|= 1.2015
At iterate 16 f = 1798.8 |proj g|= 1.4287
At iterate 17 f = 1798.3 |proj g|= 1.9001
At iterate 18 f = 1797.6 |proj g|= 1.8921
At iterate 19 f = 1796.9 |proj g|= 1.8927
At iterate 20 f = 1796.5 |proj g|= 0.9502
At iterate 21 f = 1795.9 |proj g|= 1.2526
At iterate 22 f = 1795.6 |proj g|= 1.8954
At iterate 23 f = 1795.4 |proj g|= 1.5256
At iterate 24 f = 1794.9 |proj g|= 1.8901
At iterate 25 f = 1794.7 |proj g|= 1.5025
At iterate 26 f = 1793.9 |proj g|= 1.3057
At iterate 27 f = 1793.7 |proj g|= 1.3204
At iterate 28 f = 1793.5 |proj g|= 1.4504
At iterate 29 f = 1793.4 |proj g|= 1.2978
At iterate 30 f = 1793.3 |proj g|= 1.2918
At iterate 31 f = 1793.1 |proj g|= 1.8873
At iterate 32 f = 1792.8 |proj g|= 1.8811
At iterate 33 f = 1792.6 |proj g|= 1.8802
At iterate 34 f = 1792.2 |proj g|= 1.7755
At iterate 35 f = 1791.9 |proj g|= 1.7764
At iterate 36 f = 1790.9 |proj g|= 1.7796
At iterate 37 f = 1790.8 |proj g|= 1.8808
At iterate 38 f = 1790.7 |proj g|= 1.2525
At iterate 39 f = 1790.7 |proj g|= 0.78214
At iterate 40 f = 1790.4 |proj g|= 0.82047
At iterate 41 f = 1790.3 |proj g|= 1.0794
At iterate 42 f = 1790.2 |proj g|= 1.3994
At iterate 43 f = 1790.1 |proj g|= 1.8787
At iterate 44 f = 1789.9 |proj g|= 1.2268
At iterate 45 f = 1789.9 |proj g|= 0.90265
At iterate 46 f = 1789.8 |proj g|= 1.8777
At iterate 47 f = 1789.8 |proj g|= 1.8783
At iterate 48 f = 1789.7 |proj g|= 1.5976
At iterate 49 f = 1789.7 |proj g|= 0.96825
At iterate 50 f = 1789.7 |proj g|= 1.7855
At iterate 51 f = 1789.5 |proj g|= 1.7958
At iterate 52 f = 1789.4 |proj g|= 1.7966
At iterate 53 f = 1789.4 |proj g|= 1.7892
At iterate 54 f = 1789.2 |proj g|= 1.7826
At iterate 55 f = 1789.2 |proj g|= 1.194
At iterate 56 f = 1789.2 |proj g|= 1.1932
At iterate 57 f = 1789.1 |proj g|= 1.1212
At iterate 58 f = 1789 |proj g|= 1.1135
At iterate 59 f = 1789 |proj g|= 1.875
At iterate 60 f = 1789 |proj g|= 1.7381
At iterate 61 f = 1788.9 |proj g|= 1.4756
At iterate 62 f = 1788.9 |proj g|= 1.7829
At iterate 63 f = 1788.8 |proj g|= 1.7807
At iterate 64 f = 1788.8 |proj g|= 1.7842
At iterate 65 f = 1788.8 |proj g|= 0.55799
At iterate 66 f = 1788.8 |proj g|= 0.54893
At iterate 67 f = 1788.8 |proj g|= 1.162
At iterate 68 f = 1788.8 |proj g|= 1.8723
At iterate 69 f = 1788.8 |proj g|= 1.0113
At iterate 70 f = 1788.8 |proj g|= 0.42081
At iterate 71 f = 1788.8 |proj g|= 0.40831
At iterate 72 f = 1788.8 |proj g|= 0.37755
At iterate 73 f = 1788.7 |proj g|= 0.55045
At iterate 74 f = 1788.7 |proj g|= 1.8706
At iterate 75 f = 1788.7 |proj g|= 0.62034
At iterate 76 f = 1788.7 |proj g|= 0.32525
At iterate 77 f = 1788.7 |proj g|= 0.30756
At iterate 78 f = 1788.7 |proj g|= 0.39267
At iterate 79 f = 1788.7 |proj g|= 0.63073
At iterate 80 f = 1788.7 |proj g|= 0.54849
At iterate 81 f = 1788.7 |proj g|= 0.37813
At iterate 82 f = 1788.7 |proj g|= 0.28969
At iterate 83 f = 1788.7 |proj g|= 0.40478
At iterate 84 f = 1788.7 |proj g|= 0.5533
At iterate 85 f = 1788.7 |proj g|= 0.68968
At iterate 86 f = 1788.7 |proj g|= 1.2444
At iterate 87 f = 1788.7 |proj g|= 0.099835
At iterate 88 f = 1788.7 |proj g|= 0.099618
At iterate 89 f = 1788.7 |proj g|= 0.1305
At iterate 90 f = 1788.7 |proj g|= 0.13061
At iterate 91 f = 1788.7 |proj g|= 0.47715
At iterate 92 f = 1788.7 |proj g|= 0.21897
At iterate 93 f = 1788.7 |proj g|= 0.1916
At iterate 94 f = 1788.7 |proj g|= 0.22821
At iterate 95 f = 1788.7 |proj g|= 1.0767
At iterate 96 f = 1788.7 |proj g|= 0.3527
At iterate 97 f = 1788.7 |proj g|= 0.33189
At iterate 98 f = 1788.7 |proj g|= 0.21486
At iterate 99 f = 1788.7 |proj g|= 0.1876
At iterate 100 f = 1788.7 |proj g|= 0.22777
At iterate 101 f = 1788.7 |proj g|= 0.53975
final value 1788.683283
stopped after 101 iterations
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1844.2 |proj g|= 1.901
At iterate 1 f = 1829.2 |proj g|= 1.7095
At iterate 2 f = 1826.3 |proj g|= 1.7222
At iterate 3 f = 1817.4 |proj g|= 1.8951
At iterate 4 f = 1810.3 |proj g|= 1.9269
At iterate 5 f = 1809.7 |proj g|= 1.3481
At iterate 6 f = 1808.7 |proj g|= 1.3548
At iterate 7 f = 1807.2 |proj g|= 1.9013
At iterate 8 f = 1805.9 |proj g|= 1.9149
At iterate 9 f = 1805.4 |proj g|= 1.8995
At iterate 10 f = 1803.8 |proj g|= 1.1077
At iterate 11 f = 1803.5 |proj g|= 1.9043
At iterate 12 f = 1803.4 |proj g|= 1.0665
At iterate 13 f = 1802.9 |proj g|= 1.9108
At iterate 14 f = 1802.7 |proj g|= 1.0578
At iterate 15 f = 1801.1 |proj g|= 1.9175
At iterate 16 f = 1800.8 |proj g|= 1.6693
At iterate 17 f = 1800.1 |proj g|= 1.6028
At iterate 18 f = 1799.4 |proj g|= 1.9058
At iterate 19 f = 1798 |proj g|= 1.8989
At iterate 20 f = 1797.2 |proj g|= 0.94213
At iterate 21 f = 1796.9 |proj g|= 1.8982
At iterate 22 f = 1796.3 |proj g|= 1.3185
At iterate 23 f = 1794.8 |proj g|= 1.524
At iterate 24 f = 1794.4 |proj g|= 1.5174
At iterate 25 f = 1793.8 |proj g|= 1.8934
At iterate 26 f = 1793.7 |proj g|= 1.8903
At iterate 27 f = 1793.6 |proj g|= 1.4863
At iterate 28 f = 1793.4 |proj g|= 1.4153
At iterate 29 f = 1793.3 |proj g|= 1.6275
At iterate 30 f = 1793.1 |proj g|= 1.7603
At iterate 31 f = 1792.9 |proj g|= 1.8829
At iterate 32 f = 1792.7 |proj g|= 1.8851
At iterate 33 f = 1792.5 |proj g|= 1.7591
At iterate 34 f = 1792.4 |proj g|= 0.74821
At iterate 35 f = 1792 |proj g|= 0.67033
At iterate 36 f = 1791.2 |proj g|= 0.88902
At iterate 37 f = 1791 |proj g|= 1.8859
At iterate 38 f = 1790.9 |proj g|= 0.8918
At iterate 39 f = 1790.9 |proj g|= 1.885
At iterate 40 f = 1790.9 |proj g|= 0.89031
At iterate 41 f = 1790.8 |proj g|= 0.88962
At iterate 42 f = 1790.8 |proj g|= 1.884
At iterate 43 f = 1790.7 |proj g|= 1.6486
At iterate 44 f = 1790.7 |proj g|= 0.87749
At iterate 45 f = 1790.5 |proj g|= 0.87374
At iterate 46 f = 1790.2 |proj g|= 0.92403
At iterate 47 f = 1790 |proj g|= 0.89554
At iterate 48 f = 1789.9 |proj g|= 0.85871
At iterate 49 f = 1789.8 |proj g|= 0.81672
At iterate 50 f = 1789.8 |proj g|= 1.7848
At iterate 51 f = 1789.8 |proj g|= 1.7893
At iterate 52 f = 1789.7 |proj g|= 1.7884
At iterate 53 f = 1789.7 |proj g|= 0.75862
At iterate 54 f = 1789.7 |proj g|= 0.7675
At iterate 55 f = 1789.6 |proj g|= 0.77272
At iterate 56 f = 1789.5 |proj g|= 1.8809
At iterate 57 f = 1789.4 |proj g|= 0.51739
At iterate 58 f = 1789.3 |proj g|= 1.8875
At iterate 59 f = 1789.2 |proj g|= 1.8826
At iterate 60 f = 1789.2 |proj g|= 1.7779
At iterate 61 f = 1789.2 |proj g|= 1.3471
At iterate 62 f = 1789.1 |proj g|= 0.54392
At iterate 63 f = 1789.1 |proj g|= 0.86088
At iterate 64 f = 1789.1 |proj g|= 1.8811
At iterate 65 f = 1789.1 |proj g|= 1.8811
At iterate 66 f = 1789 |proj g|= 1.8798
At iterate 67 f = 1788.9 |proj g|= 1.8743
At iterate 68 f = 1788.9 |proj g|= 0.80123
At iterate 69 f = 1788.9 |proj g|= 0.78927
At iterate 70 f = 1788.8 |proj g|= 1.8759
At iterate 71 f = 1788.8 |proj g|= 0.71153
At iterate 72 f = 1788.7 |proj g|= 0.52376
At iterate 73 f = 1788.7 |proj g|= 0.37546
At iterate 74 f = 1788.7 |proj g|= 0.32805
At iterate 75 f = 1788.7 |proj g|= 0.31573
At iterate 76 f = 1788.7 |proj g|= 0.29298
At iterate 77 f = 1788.7 |proj g|= 1.8724
At iterate 78 f = 1788.7 |proj g|= 0.77072
At iterate 79 f = 1788.7 |proj g|= 0.42465
At iterate 80 f = 1788.7 |proj g|= 0.30258
At iterate 81 f = 1788.7 |proj g|= 0.55782
At iterate 82 f = 1788.7 |proj g|= 0.59316
At iterate 83 f = 1788.7 |proj g|= 0.61489
At iterate 84 f = 1788.7 |proj g|= 0.50715
At iterate 85 f = 1788.7 |proj g|= 0.16702
At iterate 86 f = 1788.7 |proj g|= 0.28302
At iterate 87 f = 1788.7 |proj g|= 0.33463
At iterate 88 f = 1788.7 |proj g|= 0.67044
At iterate 89 f = 1788.7 |proj g|= 0.66924
At iterate 90 f = 1788.7 |proj g|= 0.21373
At iterate 91 f = 1788.7 |proj g|= 0.21343
At iterate 92 f = 1788.7 |proj g|= 0.12906
At iterate 93 f = 1788.7 |proj g|= 0.23661
At iterate 94 f = 1788.7 |proj g|= 0.15096
At iterate 95 f = 1788.7 |proj g|= 0.074849
At iterate 96 f = 1788.7 |proj g|= 0.089774
At iterate 97 f = 1788.7 |proj g|= 0.15525
At iterate 98 f = 1788.7 |proj g|= 0.20395
At iterate 99 f = 1788.7 |proj g|= 0.19495
At iterate 100 f = 1788.7 |proj g|= 0.12942
At iterate 101 f = 1788.7 |proj g|= 0.12935
final value 1788.682106
stopped after 101 iterations
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1844.5 |proj g|= 1.8965
At iterate 1 f = 1837.2 |proj g|= 1.4588
At iterate 2 f = 1830.9 |proj g|= 1.6454
ys=-1.833e-01 -gs= 5.897e+00, BFGS update SKIPPED
At iterate 3 f = 1824.7 |proj g|= 1.8975
At iterate 4 f = 1821.1 |proj g|= 1.8605
ys=-6.464e-01 -gs= 3.354e+00, BFGS update SKIPPED
At iterate 5 f = 1815 |proj g|= 1.9047
At iterate 6 f = 1810.3 |proj g|= 1.8804
At iterate 7 f = 1809.5 |proj g|= 1.5698
At iterate 8 f = 1808.6 |proj g|= 1.8781
At iterate 9 f = 1808 |proj g|= 1.8382
At iterate 10 f = 1807.5 |proj g|= 1.657
At iterate 11 f = 1806.9 |proj g|= 1.8304
At iterate 12 f = 1806.3 |proj g|= 1.8426
At iterate 13 f = 1804.9 |proj g|= 1.8294
At iterate 14 f = 1804.6 |proj g|= 1.6177
At iterate 15 f = 1804.3 |proj g|= 1.5993
At iterate 16 f = 1804.1 |proj g|= 1.853
At iterate 17 f = 1803.7 |proj g|= 1.568
At iterate 18 f = 1803.5 |proj g|= 1.8548
At iterate 19 f = 1798.5 |proj g|= 1.8698
At iterate 20 f = 1798 |proj g|= 1.5165
At iterate 21 f = 1796.5 |proj g|= 1.8676
At iterate 22 f = 1796.4 |proj g|= 1.5568
At iterate 23 f = 1796.2 |proj g|= 1.5642
At iterate 24 f = 1793.6 |proj g|= 1.5362
At iterate 25 f = 1793.3 |proj g|= 1.8045
At iterate 26 f = 1793.1 |proj g|= 1.8804
At iterate 27 f = 1792.9 |proj g|= 1.8841
At iterate 28 f = 1792.6 |proj g|= 1.8834
At iterate 29 f = 1792 |proj g|= 1.4467
At iterate 30 f = 1791.8 |proj g|= 1.2888
At iterate 31 f = 1790.9 |proj g|= 1.8787
At iterate 32 f = 1790.7 |proj g|= 1.0173
At iterate 33 f = 1790.5 |proj g|= 0.98947
At iterate 34 f = 1790.1 |proj g|= 1.8805
At iterate 35 f = 1790 |proj g|= 1.4888
At iterate 36 f = 1790 |proj g|= 0.70684
At iterate 37 f = 1790 |proj g|= 0.91985
At iterate 38 f = 1790 |proj g|= 1.6393
At iterate 39 f = 1789.8 |proj g|= 1.882
At iterate 40 f = 1789.7 |proj g|= 1.881
At iterate 41 f = 1789.6 |proj g|= 1.8776
At iterate 42 f = 1789.6 |proj g|= 0.76916
At iterate 43 f = 1789.5 |proj g|= 0.71971
At iterate 44 f = 1789.5 |proj g|= 1.1675
At iterate 45 f = 1789.5 |proj g|= 1.6044
At iterate 46 f = 1789.4 |proj g|= 1.8782
At iterate 47 f = 1789.2 |proj g|= 1.8816
At iterate 48 f = 1789.2 |proj g|= 1.801
At iterate 49 f = 1789.1 |proj g|= 1.7905
At iterate 50 f = 1789.1 |proj g|= 0.92407
At iterate 51 f = 1789.1 |proj g|= 0.83788
At iterate 52 f = 1789 |proj g|= 0.77994
At iterate 53 f = 1789 |proj g|= 0.60603
At iterate 54 f = 1788.9 |proj g|= 0.5705
At iterate 55 f = 1788.9 |proj g|= 1.8185
At iterate 56 f = 1788.9 |proj g|= 1.0668
At iterate 57 f = 1788.9 |proj g|= 1.0849
At iterate 58 f = 1788.9 |proj g|= 0.99779
At iterate 59 f = 1788.8 |proj g|= 0.30721
At iterate 60 f = 1788.8 |proj g|= 0.24539
At iterate 61 f = 1788.8 |proj g|= 1.0934
At iterate 62 f = 1788.8 |proj g|= 1.1857
At iterate 63 f = 1788.8 |proj g|= 0.40884
At iterate 64 f = 1788.8 |proj g|= 0.35038
At iterate 65 f = 1788.8 |proj g|= 0.57915
At iterate 66 f = 1788.8 |proj g|= 0.55331
At iterate 67 f = 1788.8 |proj g|= 0.75052
At iterate 68 f = 1788.7 |proj g|= 0.99308
At iterate 69 f = 1788.7 |proj g|= 1.7903
At iterate 70 f = 1788.7 |proj g|= 0.57204
At iterate 71 f = 1788.7 |proj g|= 0.30753
At iterate 72 f = 1788.7 |proj g|= 0.40401
At iterate 73 f = 1788.7 |proj g|= 1.4657
At iterate 74 f = 1788.7 |proj g|= 0.55761
At iterate 75 f = 1788.7 |proj g|= 0.49068
At iterate 76 f = 1788.7 |proj g|= 0.42876
At iterate 77 f = 1788.7 |proj g|= 0.46333
At iterate 78 f = 1788.7 |proj g|= 1.208
At iterate 79 f = 1788.7 |proj g|= 0.6762
At iterate 80 f = 1788.7 |proj g|= 0.49907
At iterate 81 f = 1788.7 |proj g|= 0.33404
At iterate 82 f = 1788.7 |proj g|= 0.16195
At iterate 83 f = 1788.7 |proj g|= 0.40693
At iterate 84 f = 1788.7 |proj g|= 0.19521
At iterate 85 f = 1788.7 |proj g|= 0.16413
At iterate 86 f = 1788.7 |proj g|= 0.13068
At iterate 87 f = 1788.7 |proj g|= 0.061973
At iterate 88 f = 1788.7 |proj g|= 0.1542
At iterate 89 f = 1788.7 |proj g|= 0.18834
At iterate 90 f = 1788.7 |proj g|= 0.11761
At iterate 91 f = 1788.7 |proj g|= 0.11372
At iterate 92 f = 1788.7 |proj g|= 0.13048
At iterate 93 f = 1788.7 |proj g|= 0.41079
At iterate 94 f = 1788.7 |proj g|= 0.13047
At iterate 95 f = 1788.7 |proj g|= 0.060687
At iterate 96 f = 1788.7 |proj g|= 0.061571
At iterate 97 f = 1788.7 |proj g|= 0.13531
At iterate 98 f = 1788.7 |proj g|= 0.20686
At iterate 99 f = 1788.7 |proj g|= 0.13803
At iterate 100 f = 1788.7 |proj g|= 0.44702
At iterate 101 f = 1788.7 |proj g|= 0.17965
final value 1788.680770
stopped after 101 iterations
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1844.5 |proj g|= 1.9627
At iterate 1 f = 1841.8 |proj g|= 1.5053
At iterate 2 f = 1835.2 |proj g|= 1.5189
At iterate 3 f = 1832.8 |proj g|= 1.9281
At iterate 4 f = 1828 |proj g|= 1.2229
At iterate 5 f = 1823.6 |proj g|= 1.0775
At iterate 6 f = 1822 |proj g|= 1.9408
At iterate 7 f = 1809.7 |proj g|= 1.8918
At iterate 8 f = 1803 |proj g|= 1.845
At iterate 9 f = 1802.6 |proj g|= 1.7398
At iterate 10 f = 1801.4 |proj g|= 1.2403
At iterate 11 f = 1799.9 |proj g|= 1.8317
At iterate 12 f = 1799 |proj g|= 1.8231
At iterate 13 f = 1797.6 |proj g|= 0.99073
At iterate 14 f = 1795.2 |proj g|= 1.0283
At iterate 15 f = 1793.4 |proj g|= 1.8782
At iterate 16 f = 1792.8 |proj g|= 1.7883
At iterate 17 f = 1792.7 |proj g|= 1.544
At iterate 18 f = 1792.6 |proj g|= 1.8611
At iterate 19 f = 1792.5 |proj g|= 1.8626
At iterate 20 f = 1791.9 |proj g|= 1.867
At iterate 21 f = 1791.7 |proj g|= 1.4722
At iterate 22 f = 1791.6 |proj g|= 1.456
At iterate 23 f = 1791.3 |proj g|= 1.7959
At iterate 24 f = 1791 |proj g|= 1.7963
At iterate 25 f = 1790.9 |proj g|= 1.8665
At iterate 26 f = 1790.7 |proj g|= 1.8672
At iterate 27 f = 1790.7 |proj g|= 0.78433
At iterate 28 f = 1790.6 |proj g|= 0.78289
At iterate 29 f = 1790.6 |proj g|= 0.7744
At iterate 30 f = 1790.4 |proj g|= 1.8684
At iterate 31 f = 1790.4 |proj g|= 1.7928
At iterate 32 f = 1790.3 |proj g|= 0.80236
At iterate 33 f = 1790.3 |proj g|= 0.82648
At iterate 34 f = 1790.1 |proj g|= 0.87867
At iterate 35 f = 1790 |proj g|= 0.8733
At iterate 36 f = 1789.8 |proj g|= 1.7905
At iterate 37 f = 1789.7 |proj g|= 1.8743
At iterate 38 f = 1789.7 |proj g|= 1.8715
At iterate 39 f = 1789.6 |proj g|= 0.80037
At iterate 40 f = 1789.6 |proj g|= 0.80415
At iterate 41 f = 1789.5 |proj g|= 0.80189
At iterate 42 f = 1789.5 |proj g|= 0.78967
At iterate 43 f = 1789.4 |proj g|= 0.75391
At iterate 44 f = 1789.3 |proj g|= 1.872
At iterate 45 f = 1789.3 |proj g|= 0.68717
At iterate 46 f = 1789.1 |proj g|= 1.4814
At iterate 47 f = 1789.1 |proj g|= 1.487
At iterate 48 f = 1789 |proj g|= 1.8742
At iterate 49 f = 1789 |proj g|= 1.4876
At iterate 50 f = 1788.9 |proj g|= 1.7243
At iterate 51 f = 1788.9 |proj g|= 0.49094
At iterate 52 f = 1788.9 |proj g|= 0.5966
At iterate 53 f = 1788.9 |proj g|= 0.59535
At iterate 54 f = 1788.9 |proj g|= 1.8733
At iterate 55 f = 1788.8 |proj g|= 0.93604
At iterate 56 f = 1788.8 |proj g|= 0.85211
At iterate 57 f = 1788.8 |proj g|= 1.4377
At iterate 58 f = 1788.8 |proj g|= 1.4772
At iterate 59 f = 1788.8 |proj g|= 0.78361
At iterate 60 f = 1788.8 |proj g|= 1.8383
At iterate 61 f = 1788.8 |proj g|= 0.3324
At iterate 62 f = 1788.8 |proj g|= 0.29873
At iterate 63 f = 1788.8 |proj g|= 0.94247
At iterate 64 f = 1788.8 |proj g|= 0.30321
At iterate 65 f = 1788.8 |proj g|= 0.32765
At iterate 66 f = 1788.7 |proj g|= 0.48866
At iterate 67 f = 1788.7 |proj g|= 1.1055
At iterate 68 f = 1788.7 |proj g|= 0.71111
At iterate 69 f = 1788.7 |proj g|= 0.50294
At iterate 70 f = 1788.7 |proj g|= 0.95867
At iterate 71 f = 1788.7 |proj g|= 0.49225
At iterate 72 f = 1788.7 |proj g|= 0.33465
At iterate 73 f = 1788.7 |proj g|= 0.3246
At iterate 74 f = 1788.7 |proj g|= 0.27462
At iterate 75 f = 1788.7 |proj g|= 0.22952
At iterate 76 f = 1788.7 |proj g|= 0.47635
At iterate 77 f = 1788.7 |proj g|= 0.36506
At iterate 78 f = 1788.7 |proj g|= 0.2705
At iterate 79 f = 1788.7 |proj g|= 0.35958
At iterate 80 f = 1788.7 |proj g|= 0.21413
At iterate 81 f = 1788.7 |proj g|= 0.083278
At iterate 82 f = 1788.7 |proj g|= 0.069125
At iterate 83 f = 1788.7 |proj g|= 0.044831
At iterate 84 f = 1788.7 |proj g|= 0.13009
At iterate 85 f = 1788.7 |proj g|= 0.083972
At iterate 86 f = 1788.7 |proj g|= 0.038698
At iterate 87 f = 1788.7 |proj g|= 0.055152
At iterate 88 f = 1788.7 |proj g|= 0.02741
At iterate 89 f = 1788.7 |proj g|= 0.026779
At iterate 90 f = 1788.7 |proj g|= 0.0423
At iterate 91 f = 1788.7 |proj g|= 0.065137
At iterate 92 f = 1788.7 |proj g|= 0.046154
At iterate 93 f = 1788.7 |proj g|= 0.041203
At iterate 94 f = 1788.7 |proj g|= 0.058948
At iterate 95 f = 1788.7 |proj g|= 0.015991
At iterate 96 f = 1788.7 |proj g|= 0.021177
iterations 96
function evaluations 123
segments explored during Cauchy searches 101
BFGS updates skipped 0
active bounds at final generalized Cauchy point 23
norm of the final projected gradient 0.0211773
final function value 1788.68
F = 1788.68
final value 1788.679245
converged
* The 4 best values (multistart) obtained are:
1788.683 1788.682 1788.681 1788.679
* The model corresponding to the best one (1788.679) is stored.
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 2 1.995574 2 1.998154 1.790352 2 2 2 1.998084 2 2 2 2 2 1.999201 2 1.996693 1.998361 2 2 1.995738 2 1.993455 1.995643 1.998832 2 2 2 2
- best initial criterion value(s) : -6193.86 -6193.966 -6195.06 -6196.14
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 6193.9 |proj g|= 1.9414
At iterate 1 f = 6099.8 |proj g|= 1.5695
At iterate 2 f = 6091.5 |proj g|= 1.903
At iterate 3 f = 6089.3 |proj g|= 1.8889
At iterate 4 f = 6083.7 |proj g|= 1.8194
At iterate 5 f = 6082.8 |proj g|= 1.8071
At iterate 6 f = 6079.7 |proj g|= 1.9016
At iterate 7 f = 6077.9 |proj g|= 1.8974
At iterate 8 f = 6068.1 |proj g|= 1.4921
At iterate 9 f = 6066.4 |proj g|= 1.4399
At iterate 10 f = 6064.5 |proj g|= 1.9119
At iterate 11 f = 6063.5 |proj g|= 1.5745
At iterate 12 f = 6060.4 |proj g|= 1.5219
At iterate 13 f = 6060 |proj g|= 1.5145
At iterate 14 f = 6059 |proj g|= 1.9074
At iterate 15 f = 6057.9 |proj g|= 1.8963
At iterate 16 f = 6056.2 |proj g|= 1.2182
At iterate 17 f = 6054.7 |proj g|= 1.8811
At iterate 18 f = 6053.9 |proj g|= 1.1028
At iterate 19 f = 6053.2 |proj g|= 1.0823
At iterate 20 f = 6053.1 |proj g|= 1.888
At iterate 21 f = 6052.9 |proj g|= 1.0711
At iterate 22 f = 6052.7 |proj g|= 1.4231
At iterate 23 f = 6052.5 |proj g|= 1.8888
At iterate 24 f = 6052.4 |proj g|= 1.3347
At iterate 25 f = 6052.4 |proj g|= 1.0392
At iterate 26 f = 6052.4 |proj g|= 1.8869
At iterate 27 f = 6052.3 |proj g|= 1.0399
At iterate 28 f = 6052.3 |proj g|= 1.0381
At iterate 29 f = 6052.3 |proj g|= 1.023
At iterate 30 f = 6052.1 |proj g|= 0.80119
At iterate 31 f = 6051.9 |proj g|= 0.74419
At iterate 32 f = 6051.9 |proj g|= 0.69194
At iterate 33 f = 6051.8 |proj g|= 1.7041
At iterate 34 f = 6051.8 |proj g|= 0.83632
At iterate 35 f = 6051.8 |proj g|= 0.83293
At iterate 36 f = 6051.7 |proj g|= 0.64953
At iterate 37 f = 6051.7 |proj g|= 0.48825
At iterate 38 f = 6051.7 |proj g|= 0.28695
At iterate 39 f = 6051.7 |proj g|= 1.6661
At iterate 40 f = 6051.7 |proj g|= 0.46224
At iterate 41 f = 6051.7 |proj g|= 0.15696
At iterate 42 f = 6051.7 |proj g|= 0.14672
At iterate 43 f = 6051.7 |proj g|= 0.11677
At iterate 44 f = 6051.7 |proj g|= 0.048306
At iterate 45 f = 6051.7 |proj g|= 0.048253
At iterate 46 f = 6051.7 |proj g|= 0.09721
At iterate 47 f = 6051.7 |proj g|= 0.11672
At iterate 48 f = 6051.7 |proj g|= 0.103
At iterate 49 f = 6051.7 |proj g|= 0.10676
At iterate 50 f = 6051.7 |proj g|= 0.012401
iterations 50
function evaluations 56
segments explored during Cauchy searches 62
BFGS updates skipped 0
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0124006
final function value 6051.72
F = 6051.72
final value 6051.720344
converged
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 6194 |proj g|= 1.896
At iterate 1 f = 6096.2 |proj g|= 1.9139
At iterate 2 f = 6095.7 |proj g|= 1.9117
At iterate 3 f = 6092.7 |proj g|= 1.9101
At iterate 4 f = 6092.3 |proj g|= 1.9112
At iterate 5 f = 6091.6 |proj g|= 1.9106
ys=-1.660e-03 -gs= 7.046e-01, BFGS update SKIPPED
At iterate 6 f = 6088.6 |proj g|= 1.9268
At iterate 7 f = 6080 |proj g|= 1.9345
ys=-1.615e-02 -gs= 8.492e+00, BFGS update SKIPPED
At iterate 8 f = 6077 |proj g|= 1.9229
At iterate 9 f = 6069.1 |proj g|= 1.8154
At iterate 10 f = 6068.5 |proj g|= 1.6384
At iterate 11 f = 6068 |proj g|= 1.9184
At iterate 12 f = 6067.9 |proj g|= 1.9211
At iterate 13 f = 6067 |proj g|= 1.9248
At iterate 14 f = 6063.9 |proj g|= 1.9201
At iterate 15 f = 6062.4 |proj g|= 1.9124
At iterate 16 f = 6061.1 |proj g|= 1.7501
At iterate 17 f = 6059.5 |proj g|= 1.9052
At iterate 18 f = 6059.1 |proj g|= 1.7219
At iterate 19 f = 6057.7 |proj g|= 1.7267
At iterate 20 f = 6056.9 |proj g|= 1.7136
At iterate 21 f = 6056.1 |proj g|= 1.6472
At iterate 22 f = 6055.3 |proj g|= 1.8976
At iterate 23 f = 6055.2 |proj g|= 1.6393
At iterate 24 f = 6055 |proj g|= 1.629
At iterate 25 f = 6054.6 |proj g|= 1.8972
At iterate 26 f = 6054.6 |proj g|= 1.5981
At iterate 27 f = 6054.5 |proj g|= 1.5928
At iterate 28 f = 6054.4 |proj g|= 1.5774
At iterate 29 f = 6054.2 |proj g|= 1.8957
At iterate 30 f = 6054.1 |proj g|= 1.8949
At iterate 31 f = 6052.6 |proj g|= 1.2987
At iterate 32 f = 6052.3 |proj g|= 1.2399
At iterate 33 f = 6051.9 |proj g|= 0.96675
At iterate 34 f = 6051.8 |proj g|= 0.96791
At iterate 35 f = 6051.8 |proj g|= 0.94212
At iterate 36 f = 6051.8 |proj g|= 1.8852
At iterate 37 f = 6051.8 |proj g|= 0.71337
At iterate 38 f = 6051.8 |proj g|= 0.71134
At iterate 39 f = 6051.7 |proj g|= 0.3276
At iterate 40 f = 6051.7 |proj g|= 0.19604
At iterate 41 f = 6051.7 |proj g|= 0.14445
At iterate 42 f = 6051.7 |proj g|= 0.13129
At iterate 43 f = 6051.7 |proj g|= 0.11073
At iterate 44 f = 6051.7 |proj g|= 0.65001
At iterate 45 f = 6051.7 |proj g|= 0.50626
At iterate 46 f = 6051.7 |proj g|= 0.21399
At iterate 47 f = 6051.7 |proj g|= 0.1184
At iterate 48 f = 6051.7 |proj g|= 0.11665
At iterate 49 f = 6051.7 |proj g|= 0.076832
At iterate 50 f = 6051.7 |proj g|= 0.11684
At iterate 51 f = 6051.7 |proj g|= 0.1167
At iterate 52 f = 6051.7 |proj g|= 0.09314
At iterate 53 f = 6051.7 |proj g|= 0.10182
At iterate 54 f = 6051.7 |proj g|= 0.050777
At iterate 55 f = 6051.7 |proj g|= 0.025595
At iterate 56 f = 6051.7 |proj g|= 0.11416
iterations 56
function evaluations 73
segments explored during Cauchy searches 68
BFGS updates skipped 2
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.114159
final function value 6051.72
F = 6051.72
final value 6051.720361
converged
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 6195.1 |proj g|= 1.8444
At iterate 1 f = 6190.4 |proj g|= 1.567
At iterate 2 f = 6180.2 |proj g|= 1.8736
At iterate 3 f = 6172.8 |proj g|= 1.8258
ys=-8.212e-01 -gs= 6.952e+00, BFGS update SKIPPED
At iterate 4 f = 6153.1 |proj g|= 1.7389
ys=-2.505e+00 -gs= 1.744e+01, BFGS update SKIPPED
At iterate 5 f = 6141.1 |proj g|= 1.694
At iterate 6 f = 6130.4 |proj g|= 1.9376
At iterate 7 f = 6113 |proj g|= 1.8015
At iterate 8 f = 6107.7 |proj g|= 1.8026
At iterate 9 f = 6085.2 |proj g|= 1.6781
At iterate 10 f = 6084 |proj g|= 1.8835
At iterate 11 f = 6082.6 |proj g|= 1.8806
At iterate 12 f = 6072.2 |proj g|= 1.8738
At iterate 13 f = 6070.1 |proj g|= 1.4886
At iterate 14 f = 6065.8 |proj g|= 1.381
At iterate 15 f = 6064.4 |proj g|= 1.3835
At iterate 16 f = 6064 |proj g|= 1.8899
At iterate 17 f = 6062.2 |proj g|= 1.8814
At iterate 18 f = 6061.5 |proj g|= 1.421
At iterate 19 f = 6059.1 |proj g|= 1.4174
At iterate 20 f = 6058.3 |proj g|= 1.4164
At iterate 21 f = 6056.7 |proj g|= 1.4007
At iterate 22 f = 6055.7 |proj g|= 1.901
At iterate 23 f = 6055.3 |proj g|= 1.4749
At iterate 24 f = 6054.9 |proj g|= 1.4674
At iterate 25 f = 6054.7 |proj g|= 1.8948
At iterate 26 f = 6054.5 |proj g|= 1.8955
At iterate 27 f = 6054.4 |proj g|= 1.8936
At iterate 28 f = 6053.8 |proj g|= 1.1737
At iterate 29 f = 6053.3 |proj g|= 1.8945
At iterate 30 f = 6053.2 |proj g|= 1.8917
At iterate 31 f = 6053.1 |proj g|= 1.353
At iterate 32 f = 6052.9 |proj g|= 1.4434
At iterate 33 f = 6052.8 |proj g|= 1.8925
At iterate 34 f = 6052.7 |proj g|= 1.441
At iterate 35 f = 6052.4 |proj g|= 1.4288
At iterate 36 f = 6052.3 |proj g|= 1.0089
At iterate 37 f = 6052.2 |proj g|= 1.8889
At iterate 38 f = 6052.1 |proj g|= 1.8888
At iterate 39 f = 6052.1 |proj g|= 1.6079
At iterate 40 f = 6052 |proj g|= 1.2944
At iterate 41 f = 6052 |proj g|= 1.1205
At iterate 42 f = 6052 |proj g|= 1.8863
At iterate 43 f = 6052 |proj g|= 1.429
At iterate 44 f = 6051.9 |proj g|= 0.92281
At iterate 45 f = 6051.9 |proj g|= 1.3317
At iterate 46 f = 6051.9 |proj g|= 1.3044
At iterate 47 f = 6051.9 |proj g|= 0.78826
At iterate 48 f = 6051.9 |proj g|= 0.80595
At iterate 49 f = 6051.9 |proj g|= 0.91726
At iterate 50 f = 6051.8 |proj g|= 1.8835
At iterate 51 f = 6051.7 |proj g|= 1.8839
At iterate 52 f = 6051.7 |proj g|= 0.62815
At iterate 53 f = 6051.7 |proj g|= 0.59062
At iterate 54 f = 6051.7 |proj g|= 0.09207
At iterate 55 f = 6051.7 |proj g|= 0.069011
At iterate 56 f = 6051.7 |proj g|= 0.016776
At iterate 57 f = 6051.7 |proj g|= 0.007937
iterations 57
function evaluations 64
segments explored during Cauchy searches 65
BFGS updates skipped 2
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.00793695
final function value 6051.72
F = 6051.72
final value 6051.720335
converged
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 6196.1 |proj g|= 1.9441
At iterate 1 f = 6158.6 |proj g|= 1.8479
At iterate 2 f = 6131.1 |proj g|= 1.8288
ys=-4.125e+00 -gs= 2.537e+01, BFGS update SKIPPED
At iterate 3 f = 6093.3 |proj g|= 1.7964
At iterate 4 f = 6091.6 |proj g|= 1.8966
At iterate 5 f = 6089.7 |proj g|= 1.7487
At iterate 6 f = 6078.4 |proj g|= 1.7556
At iterate 7 f = 6069.5 |proj g|= 1.7184
At iterate 8 f = 6066.2 |proj g|= 1.6136
At iterate 9 f = 6065.3 |proj g|= 1.9357
At iterate 10 f = 6063.3 |proj g|= 1.9207
At iterate 11 f = 6061.3 |proj g|= 1.909
At iterate 12 f = 6059.9 |proj g|= 1.8995
At iterate 13 f = 6057.7 |proj g|= 1.7114
At iterate 14 f = 6057.6 |proj g|= 1.7106
At iterate 15 f = 6057.1 |proj g|= 1.9015
At iterate 16 f = 6056.4 |proj g|= 1.9085
At iterate 17 f = 6056 |proj g|= 1.8998
At iterate 18 f = 6055.9 |proj g|= 1.6836
At iterate 19 f = 6055.7 |proj g|= 1.6795
At iterate 20 f = 6055.6 |proj g|= 1.8978
At iterate 21 f = 6055.5 |proj g|= 1.8959
At iterate 22 f = 6055.3 |proj g|= 1.6477
At iterate 23 f = 6055 |proj g|= 1.6138
At iterate 24 f = 6054.5 |proj g|= 1.5619
At iterate 25 f = 6054.1 |proj g|= 1.8977
At iterate 26 f = 6054 |proj g|= 1.8975
At iterate 27 f = 6053.4 |proj g|= 1.4395
At iterate 28 f = 6053 |proj g|= 1.889
At iterate 29 f = 6052.3 |proj g|= 1.8807
At iterate 30 f = 6052 |proj g|= 0.72993
At iterate 31 f = 6051.9 |proj g|= 0.84726
At iterate 32 f = 6051.8 |proj g|= 1.8843
At iterate 33 f = 6051.8 |proj g|= 0.47819
At iterate 34 f = 6051.8 |proj g|= 1.8846
At iterate 35 f = 6051.7 |proj g|= 0.70317
At iterate 36 f = 6051.7 |proj g|= 0.58719
At iterate 37 f = 6051.7 |proj g|= 0.53426
At iterate 38 f = 6051.7 |proj g|= 0.28181
At iterate 39 f = 6051.7 |proj g|= 0.27471
At iterate 40 f = 6051.7 |proj g|= 0.71023
At iterate 41 f = 6051.7 |proj g|= 0.7442
At iterate 42 f = 6051.7 |proj g|= 0.17579
At iterate 43 f = 6051.7 |proj g|= 0.12481
At iterate 44 f = 6051.7 |proj g|= 0.11666
At iterate 45 f = 6051.7 |proj g|= 0.80637
At iterate 46 f = 6051.7 |proj g|= 0.047588
At iterate 47 f = 6051.7 |proj g|= 0.041069
At iterate 48 f = 6051.7 |proj g|= 0.1168
At iterate 49 f = 6051.7 |proj g|= 0.11675
At iterate 50 f = 6051.7 |proj g|= 0.043049
At iterate 51 f = 6051.7 |proj g|= 0.042103
At iterate 52 f = 6051.7 |proj g|= 0.065384
At iterate 53 f = 6051.7 |proj g|= 0.11522
At iterate 54 f = 6051.7 |proj g|= 0.044015
iterations 54
function evaluations 68
segments explored during Cauchy searches 64
BFGS updates skipped 1
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0440146
final function value 6051.72
F = 6051.72
final value 6051.720388
converged
* The 4 best values (multistart) obtained are:
6051.72 6051.72 6051.72 6051.72
* The model corresponding to the best one (6051.72) is stored.
[1] 10.81038
[1] 369550.6
[1] 8.353776
[1] 5.220835
The problem with doing constraint first is that you end up with a non-hypercube shaped input space (corners have been knocked off by constraint), which is not ideal. We might therefore want different sensitivity measures for pre- and post-constrained ensemble.
[1] 8.353776
[1] 8.353615
[1] 5.220835
[1] 5.220903
It looks as though bwl_io is very influential across a number of variables, even though it didn’t appear that interesting in the parginal plots. Why is that?
A closer look at bwl_io now that the impact of f0_io has been removed shows a large number of “zero” NPP when bwl_io is at low values, which could well be the source of apparent sensitivity.
Take only higher values of bwl_io for the next round of constraints.
# Plot the regular km emulator.
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -5137.861
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 5137.9 |proj g|= 1.8356
At iterate 1 f = 5135.4 |proj g|= 1.7202
At iterate 2 f = 5131.9 |proj g|= 1.6
At iterate 3 f = 5131.2 |proj g|= 1.6826
At iterate 4 f = 5122.2 |proj g|= 1.6838
At iterate 5 f = 5119.7 |proj g|= 1.6697
At iterate 6 f = 5110.5 |proj g|= 1.6667
At iterate 7 f = 5105.6 |proj g|= 1.3939
At iterate 8 f = 5098.3 |proj g|= 1.4145
At iterate 9 f = 5093.3 |proj g|= 1.592
At iterate 10 f = 5091.8 |proj g|= 1.5626
At iterate 11 f = 5089.7 |proj g|= 1.4986
At iterate 12 f = 5088.5 |proj g|= 1.3214
At iterate 13 f = 5085.1 |proj g|= 1.501
At iterate 14 f = 5083.9 |proj g|= 1.4621
At iterate 15 f = 5082.1 |proj g|= 1.4517
At iterate 16 f = 5081.1 |proj g|= 1.494
At iterate 17 f = 5080.6 |proj g|= 0.95799
At iterate 18 f = 5080.4 |proj g|= 1.3623
At iterate 19 f = 5080 |proj g|= 0.56793
At iterate 20 f = 5079.9 |proj g|= 1.4786
At iterate 21 f = 5079.8 |proj g|= 0.46896
At iterate 22 f = 5079.8 |proj g|= 0.28689
At iterate 23 f = 5079.8 |proj g|= 0.57427
At iterate 24 f = 5079.7 |proj g|= 0.57268
At iterate 25 f = 5079.7 |proj g|= 0.36422
At iterate 26 f = 5079.7 |proj g|= 0.16654
At iterate 27 f = 5079.7 |proj g|= 0.22658
At iterate 28 f = 5079.7 |proj g|= 0.27311
At iterate 29 f = 5079.7 |proj g|= 0.21159
At iterate 30 f = 5079.7 |proj g|= 0.52762
At iterate 31 f = 5079.7 |proj g|= 0.2638
At iterate 32 f = 5079.7 |proj g|= 0.15332
At iterate 33 f = 5079.7 |proj g|= 0.16807
At iterate 34 f = 5079.7 |proj g|= 0.2156
At iterate 35 f = 5079.7 |proj g|= 0.093168
At iterate 36 f = 5079.7 |proj g|= 0.069853
At iterate 37 f = 5079.7 |proj g|= 0.11439
At iterate 38 f = 5079.7 |proj g|= 0.18854
At iterate 39 f = 5079.7 |proj g|= 0.14403
At iterate 40 f = 5079.6 |proj g|= 0.075684
At iterate 41 f = 5079.6 |proj g|= 0.045529
At iterate 42 f = 5079.6 |proj g|= 0.029163
At iterate 43 f = 5079.6 |proj g|= 0.028335
At iterate 44 f = 5079.6 |proj g|= 0.081292
At iterate 45 f = 5079.6 |proj g|= 0.03851
At iterate 46 f = 5079.6 |proj g|= 0.011842
At iterate 47 f = 5079.6 |proj g|= 0.010565
iterations 47
function evaluations 54
segments explored during Cauchy searches 60
BFGS updates skipped 0
active bounds at final generalized Cauchy point 20
norm of the final projected gradient 0.0105653
final function value 5079.65
F = 5079.65
final value 5079.648744
converged
It looks as though both the absolute value and the change over time are controlled by the same parameters.
The emulators appear to be at least capturing the broad response for all of the output variables.
First, plot the straight kriging emulators
[1] 5.433077
[1] 5.142409
[1] 5.2209
[1] 4.193005
[1] 4.219396
[1] 4.531841
[1] 5.629425
[1] 8.565018
[1] 5.224583
[1] 9.223968
[1] 7.089
[1] 7.78196
[1] 6.900971
[1] 6.929568
Next, plot the twostep glmnet/km emulators
We use the leave-one-out Mean Absolute Error, expressed as a percentage of the range of the output across the ensemble. We find that the twostep emulatorisn’t significantly more accurate, and is indeed less accurate for tree fraction.
[1] 5.410355
[1] 4.733929
[1] 5.215455
[1] 4.173501
[1] 4.03861
[1] 4.597044
[1] 5.501413
[1] 8.563912
[1] 5.213222
[1] 11.36184
[1] 7.260354
[1] 7.707121
[1] 6.725096
[1] 6.865819
To use History Matching, we need to specify targets for various model outputs. We treat these targets as “observations” with an uncertainty where a model run marked as “implausible” (beyond 3sd) matches the hard boundaries previously identified by A. Wiltshire as being desirable.
Choose the centre of the (implied) uniform distribution. cs_gb.target = (3000 - 750) / 2 = 1125 cv.target = (800 - 300) / 2 = 250 npp_n_gb.target = (80 - 35) / 2 = 22.5
nbp.target = 0 (gpp.target = 75) (runoff.target = 1)
(to do: visualise implausibility of design and loo emulated values of design as histogram)
Plot the original ensemble, plus the stuff that matches the constraints. Include anomalies.
We find that adding a constraint to level 2 constrains the behaviour of the soil carbon pool, but not the vegetation carbon pool. Probably not a surprise, given the Soil Carbon pool (and changes) are much larger than the vegetation carbon pool.
We plot the JULES total land carbon change from 1850, together with “observations” from the Global Carbon Budget (Friedligstein et al, 2020). The “observations are the cumulative sum of the ‘land sink’ column on the”historical observations tab of the data from the paper. It states “The land sink is the average of several dynamic global vegetation models that reproduce the observed mean total land sink of the 1990s.”
I think to make a fair comparison, we might need to take the diff of the anomaly (or just the timeseries), as the above represents the total change since 1850, and I think all the “land sink” stuff looks at the year-to-year changes.
What do these parts of parameter space tell us about what we need to do to have the model reproduce reality?
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -106.4158
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 106.42 |proj g|= 1.8453
At iterate 1 f = 103.09 |proj g|= 1.819
At iterate 2 f = 97.125 |proj g|= 1.6466
At iterate 3 f = 95.132 |proj g|= 1.6743
At iterate 4 f = 91.572 |proj g|= 1.6622
At iterate 5 f = 89.835 |proj g|= 1.6389
At iterate 6 f = 88.807 |proj g|= 1.5108
At iterate 7 f = 87.513 |proj g|= 1.7187
At iterate 8 f = 86.326 |proj g|= 1.638
At iterate 9 f = 85.522 |proj g|= 1.5586
At iterate 10 f = 85.21 |proj g|= 1.342
At iterate 11 f = 84.63 |proj g|= 1.2524
At iterate 12 f = 84.111 |proj g|= 1.2971
At iterate 13 f = 83.72 |proj g|= 1.1355
At iterate 14 f = 83.52 |proj g|= 1.0089
At iterate 15 f = 82.998 |proj g|= 0.69506
At iterate 16 f = 82.939 |proj g|= 1.7322
At iterate 17 f = 82.779 |proj g|= 1.7049
At iterate 18 f = 82.708 |proj g|= 0.80111
At iterate 19 f = 82.614 |proj g|= 0.71523
At iterate 20 f = 82.575 |proj g|= 0.65009
At iterate 21 f = 82.517 |proj g|= 0.61643
At iterate 22 f = 82.482 |proj g|= 0.53998
At iterate 23 f = 82.35 |proj g|= 1.4933
At iterate 24 f = 82.31 |proj g|= 1.4399
At iterate 25 f = 82.247 |proj g|= 0.83833
At iterate 26 f = 82.179 |proj g|= 0.6743
At iterate 27 f = 82.125 |proj g|= 0.61099
At iterate 28 f = 82.077 |proj g|= 0.69756
At iterate 29 f = 82.068 |proj g|= 0.50643
At iterate 30 f = 82.045 |proj g|= 0.53157
At iterate 31 f = 81.976 |proj g|= 1.2632
At iterate 32 f = 81.806 |proj g|= 1.4306
At iterate 33 f = 81.751 |proj g|= 1.1701
At iterate 34 f = 81.628 |proj g|= 1.1438
At iterate 35 f = 81.556 |proj g|= 1.6474
At iterate 36 f = 81.483 |proj g|= 1.4191
At iterate 37 f = 81.399 |proj g|= 0.67894
At iterate 38 f = 81.324 |proj g|= 0.33289
At iterate 39 f = 81.314 |proj g|= 1.0959
At iterate 40 f = 81.302 |proj g|= 0.32766
At iterate 41 f = 81.295 |proj g|= 0.3297
At iterate 42 f = 81.277 |proj g|= 0.36728
At iterate 43 f = 81.265 |proj g|= 0.30341
At iterate 44 f = 81.247 |proj g|= 0.69733
At iterate 45 f = 81.216 |proj g|= 1.1062
At iterate 46 f = 81.177 |proj g|= 1.1638
At iterate 47 f = 81.155 |proj g|= 0.71712
At iterate 48 f = 81.109 |proj g|= 0.51864
At iterate 49 f = 81.077 |proj g|= 0.41088
At iterate 50 f = 81.059 |proj g|= 0.39031
At iterate 51 f = 81.055 |proj g|= 0.33235
At iterate 52 f = 81.052 |proj g|= 0.14717
At iterate 53 f = 81.049 |proj g|= 0.1497
At iterate 54 f = 81.043 |proj g|= 0.15788
At iterate 55 f = 81.035 |proj g|= 0.29801
At iterate 56 f = 81.033 |proj g|= 0.57956
At iterate 57 f = 81.022 |proj g|= 0.39169
At iterate 58 f = 81.018 |proj g|= 0.16287
At iterate 59 f = 81.015 |proj g|= 0.16427
At iterate 60 f = 81.003 |proj g|= 0.17202
At iterate 61 f = 80.998 |proj g|= 0.55788
At iterate 62 f = 80.983 |proj g|= 0.3853
At iterate 63 f = 80.978 |proj g|= 0.14235
At iterate 64 f = 80.973 |proj g|= 0.16257
At iterate 65 f = 80.969 |proj g|= 0.55124
At iterate 66 f = 80.966 |proj g|= 0.17098
At iterate 67 f = 80.964 |proj g|= 0.15834
At iterate 68 f = 80.961 |proj g|= 0.19492
At iterate 69 f = 80.96 |proj g|= 0.085999
At iterate 70 f = 80.959 |proj g|= 0.20316
At iterate 71 f = 80.957 |proj g|= 0.22804
At iterate 72 f = 80.956 |proj g|= 0.18647
At iterate 73 f = 80.955 |proj g|= 0.058348
At iterate 74 f = 80.955 |proj g|= 0.054364
At iterate 75 f = 80.953 |proj g|= 0.091177
At iterate 76 f = 80.953 |proj g|= 0.31904
At iterate 77 f = 80.951 |proj g|= 0.21525
At iterate 78 f = 80.947 |proj g|= 0.12164
At iterate 79 f = 80.943 |proj g|= 0.12463
At iterate 80 f = 80.941 |proj g|= 0.18732
At iterate 81 f = 80.939 |proj g|= 0.08267
At iterate 82 f = 80.938 |proj g|= 0.24049
At iterate 83 f = 80.938 |proj g|= 0.17522
At iterate 84 f = 80.937 |proj g|= 0.10918
At iterate 85 f = 80.937 |proj g|= 0.079848
At iterate 86 f = 80.936 |proj g|= 0.045969
At iterate 87 f = 80.936 |proj g|= 0.043919
At iterate 88 f = 80.936 |proj g|= 0.12858
At iterate 89 f = 80.936 |proj g|= 0.04115
At iterate 90 f = 80.936 |proj g|= 0.025463
At iterate 91 f = 80.936 |proj g|= 0.018393
At iterate 92 f = 80.936 |proj g|= 0.020207
At iterate 93 f = 80.935 |proj g|= 0.050297
At iterate 94 f = 80.935 |proj g|= 0.018193
At iterate 95 f = 80.935 |proj g|= 0.015472
At iterate 96 f = 80.935 |proj g|= 0.014675
At iterate 97 f = 80.935 |proj g|= 0.072348
At iterate 98 f = 80.935 |proj g|= 0.02166
At iterate 99 f = 80.935 |proj g|= 0.017904
At iterate 100 f = 80.935 |proj g|= 0.040947
At iterate 101 f = 80.935 |proj g|= 0.024479
final value 80.934858
stopped after 101 iterations
The function addNroyDesignPoints builds an emulator for each model output in Y. It compares the output of each emulator at a number of candidate desin points, and chooses a space-filling set of them that that are Not Ruled Out Yet (statistically close to the observation at Y_target).
The function write_jules_design here simply takes the points calculated by addNroyDesignPoints and writes them to configuration files.
Check that the augmented design produces what we expect. New ensemble members should be somewhat constrained within the boundaries of the original design, if the comparison to data offers any constraint.
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -109.2909
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 109.29 |proj g|= 1.871
At iterate 1 f = 100.41 |proj g|= 1.3975
At iterate 2 f = 96.533 |proj g|= 1.3088
At iterate 3 f = 95.227 |proj g|= 1.2724
At iterate 4 f = 92.512 |proj g|= 1.5964
At iterate 5 f = 91.143 |proj g|= 1.7368
At iterate 6 f = 89.279 |proj g|= 1.6991
At iterate 7 f = 88.195 |proj g|= 1.7046
At iterate 8 f = 87.014 |proj g|= 1.5114
At iterate 9 f = 85.689 |proj g|= 1.0118
At iterate 10 f = 85.023 |proj g|= 1.7091
At iterate 11 f = 84.569 |proj g|= 1.3563
At iterate 12 f = 84.158 |proj g|= 0.80825
At iterate 13 f = 84.09 |proj g|= 1.6904
At iterate 14 f = 83.823 |proj g|= 1.6909
At iterate 15 f = 83.55 |proj g|= 0.81435
At iterate 16 f = 83.371 |proj g|= 0.96783
At iterate 17 f = 83.094 |proj g|= 1.1149
At iterate 18 f = 82.902 |proj g|= 0.63225
At iterate 19 f = 82.876 |proj g|= 1.6631
At iterate 20 f = 82.821 |proj g|= 0.47627
At iterate 21 f = 82.812 |proj g|= 0.44135
At iterate 22 f = 82.724 |proj g|= 0.87459
At iterate 23 f = 82.656 |proj g|= 0.80084
At iterate 24 f = 82.618 |proj g|= 0.37465
At iterate 25 f = 82.608 |proj g|= 0.27736
At iterate 26 f = 82.599 |proj g|= 0.18904
At iterate 27 f = 82.587 |proj g|= 0.29552
At iterate 28 f = 82.563 |proj g|= 0.53817
At iterate 29 f = 82.536 |proj g|= 0.389
At iterate 30 f = 82.52 |proj g|= 0.27542
At iterate 31 f = 82.509 |proj g|= 0.43571
At iterate 32 f = 82.496 |proj g|= 0.39767
At iterate 33 f = 82.472 |proj g|= 0.67351
At iterate 34 f = 82.447 |proj g|= 1.0541
At iterate 35 f = 82.426 |proj g|= 0.41252
At iterate 36 f = 82.407 |proj g|= 0.36104
At iterate 37 f = 82.39 |proj g|= 0.37105
At iterate 38 f = 82.352 |proj g|= 0.40153
At iterate 39 f = 82.345 |proj g|= 0.48576
At iterate 40 f = 82.338 |proj g|= 0.30354
At iterate 41 f = 82.333 |proj g|= 0.28668
At iterate 42 f = 82.321 |proj g|= 0.24732
At iterate 43 f = 82.309 |proj g|= 0.34561
At iterate 44 f = 82.295 |proj g|= 0.29314
At iterate 45 f = 82.283 |proj g|= 0.5463
At iterate 46 f = 82.269 |proj g|= 0.55403
At iterate 47 f = 82.266 |proj g|= 0.38142
At iterate 48 f = 82.264 |proj g|= 0.15815
At iterate 49 f = 82.26 |proj g|= 0.11062
At iterate 50 f = 82.258 |proj g|= 0.1582
At iterate 51 f = 82.255 |proj g|= 0.22033
At iterate 52 f = 82.254 |proj g|= 0.13105
At iterate 53 f = 82.253 |proj g|= 0.05804
At iterate 54 f = 82.253 |proj g|= 0.15167
At iterate 55 f = 82.253 |proj g|= 0.054998
At iterate 56 f = 82.252 |proj g|= 0.03075
At iterate 57 f = 82.252 |proj g|= 0.029696
At iterate 58 f = 82.252 |proj g|= 0.016543
At iterate 59 f = 82.252 |proj g|= 0.014205
At iterate 60 f = 82.252 |proj g|= 0.026848
At iterate 61 f = 82.252 |proj g|= 0.027097
At iterate 62 f = 82.252 |proj g|= 0.012538
At iterate 63 f = 82.252 |proj g|= 0.011962
At iterate 64 f = 82.252 |proj g|= 0.011229
At iterate 65 f = 82.252 |proj g|= 0.01993
At iterate 66 f = 82.252 |proj g|= 0.012827
At iterate 67 f = 82.252 |proj g|= 0.016769
At iterate 68 f = 82.252 |proj g|= 0.038307
At iterate 69 f = 82.252 |proj g|= 0.02021
At iterate 70 f = 82.252 |proj g|= 0.01395
At iterate 71 f = 82.252 |proj g|= 0.019049
At iterate 72 f = 82.252 |proj g|= 0.028594
At iterate 73 f = 82.252 |proj g|= 0.033337
At iterate 74 f = 82.252 |proj g|= 0.020535
At iterate 75 f = 82.252 |proj g|= 0.0075909
At iterate 76 f = 82.252 |proj g|= 0.0055178
At iterate 77 f = 82.252 |proj g|= 0.033389
At iterate 78 f = 82.252 |proj g|= 0.0071442
At iterate 79 f = 82.252 |proj g|= 0.0111
At iterate 80 f = 82.252 |proj g|= 0.01729
At iterate 81 f = 82.252 |proj g|= 0.01571
At iterate 82 f = 82.252 |proj g|= 0.029493
At iterate 83 f = 82.252 |proj g|= 0.0070285
At iterate 84 f = 82.252 |proj g|= 0.0084202
At iterate 85 f = 82.252 |proj g|= 0.012163
At iterate 86 f = 82.252 |proj g|= 0.089832
At iterate 87 f = 82.252 |proj g|= 0.03557
At iterate 88 f = 82.252 |proj g|= 0.010387
At iterate 89 f = 82.252 |proj g|= 0.02101
At iterate 90 f = 82.252 |proj g|= 0.017472
At iterate 91 f = 82.252 |proj g|= 0.025682
At iterate 92 f = 82.252 |proj g|= 0.0079639
At iterate 93 f = 82.252 |proj g|= 0.0073904
At iterate 94 f = 82.252 |proj g|= 0.018008
At iterate 95 f = 82.252 |proj g|= 0.0062852
At iterate 96 f = 82.252 |proj g|= 0.0046032
At iterate 97 f = 82.252 |proj g|= 0.0028965
At iterate 98 f = 82.251 |proj g|= 0.0083294
At iterate 99 f = 82.251 |proj g|= 0.004218
At iterate 100 f = 82.251 |proj g|= 0.0024721
At iterate 101 f = 82.251 |proj g|= 0.0026606
final value 82.251495
stopped after 101 iterations
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -1383.375
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1383.4 |proj g|= 1.7664
At iterate 1 f = 1372.9 |proj g|= 1.6648
At iterate 2 f = 1364.8 |proj g|= 1.7747
At iterate 3 f = 1356.7 |proj g|= 1.8961
At iterate 4 f = 1350.1 |proj g|= 1.7144
At iterate 5 f = 1347.8 |proj g|= 1.6582
At iterate 6 f = 1341.7 |proj g|= 1.481
At iterate 7 f = 1337.6 |proj g|= 1.6235
At iterate 8 f = 1335.4 |proj g|= 1.4778
At iterate 9 f = 1334 |proj g|= 0.92855
At iterate 10 f = 1332 |proj g|= 1.5095
At iterate 11 f = 1330.6 |proj g|= 1.0917
At iterate 12 f = 1329.3 |proj g|= 1.2354
At iterate 13 f = 1328.5 |proj g|= 0.82781
At iterate 14 f = 1328.3 |proj g|= 0.9746
At iterate 15 f = 1328.1 |proj g|= 0.55288
At iterate 16 f = 1328 |proj g|= 0.39538
At iterate 17 f = 1327.8 |proj g|= 0.30206
At iterate 18 f = 1327.8 |proj g|= 1.4645
At iterate 19 f = 1327.8 |proj g|= 0.56521
At iterate 20 f = 1327.8 |proj g|= 0.21692
At iterate 21 f = 1327.8 |proj g|= 0.56111
At iterate 22 f = 1327.8 |proj g|= 0.51424
At iterate 23 f = 1327.8 |proj g|= 0.53842
At iterate 24 f = 1327.8 |proj g|= 0.1763
At iterate 25 f = 1327.7 |proj g|= 0.20315
At iterate 26 f = 1327.7 |proj g|= 0.14516
At iterate 27 f = 1327.7 |proj g|= 0.054732
At iterate 28 f = 1327.7 |proj g|= 0.084036
At iterate 29 f = 1327.7 |proj g|= 0.045267
At iterate 30 f = 1327.7 |proj g|= 0.046735
At iterate 31 f = 1327.7 |proj g|= 0.14927
At iterate 32 f = 1327.7 |proj g|= 0.11387
At iterate 33 f = 1327.7 |proj g|= 0.025094
At iterate 34 f = 1327.7 |proj g|= 0.018372
At iterate 35 f = 1327.7 |proj g|= 0.018539
At iterate 36 f = 1327.7 |proj g|= 0.017334
At iterate 37 f = 1327.7 |proj g|= 0.035245
At iterate 38 f = 1327.7 |proj g|= 0.016258
At iterate 39 f = 1327.7 |proj g|= 0.0075892
At iterate 40 f = 1327.7 |proj g|= 0.008652
At iterate 41 f = 1327.7 |proj g|= 0.023457
iterations 41
function evaluations 48
segments explored during Cauchy searches 62
BFGS updates skipped 0
active bounds at final generalized Cauchy point 20
norm of the final projected gradient 0.0234567
final function value 1327.74
F = 1327.74
final value 1327.742870
converged
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -2561.221
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 2561.2 |proj g|= 1.8661
At iterate 1 f = 2533.6 |proj g|= 2
At iterate 2 f = 2519.6 |proj g|= 2
At iterate 3 f = 2516.5 |proj g|= 1.9534
At iterate 4 f = 2510.4 |proj g|= 1.7867
At iterate 5 f = 2505 |proj g|= 1.6826
At iterate 6 f = 2501.5 |proj g|= 1.618
At iterate 7 f = 2494.2 |proj g|= 1.3844
At iterate 8 f = 2492.2 |proj g|= 1.6032
At iterate 9 f = 2491.3 |proj g|= 1.589
At iterate 10 f = 2488.9 |proj g|= 1.53
At iterate 11 f = 2487.5 |proj g|= 1.0728
At iterate 12 f = 2487.2 |proj g|= 1.0711
At iterate 13 f = 2487 |proj g|= 1.1027
At iterate 14 f = 2486.8 |proj g|= 0.47609
At iterate 15 f = 2486.8 |proj g|= 0.42024
At iterate 16 f = 2486.8 |proj g|= 1.2305
At iterate 17 f = 2486.8 |proj g|= 0.38381
At iterate 18 f = 2486.8 |proj g|= 0.22702
At iterate 19 f = 2486.8 |proj g|= 0.21667
At iterate 20 f = 2486.8 |proj g|= 0.17646
At iterate 21 f = 2486.7 |proj g|= 0.202
At iterate 22 f = 2486.7 |proj g|= 0.44121
At iterate 23 f = 2486.7 |proj g|= 0.16229
At iterate 24 f = 2486.7 |proj g|= 0.066349
At iterate 25 f = 2486.7 |proj g|= 0.018863
At iterate 26 f = 2486.7 |proj g|= 0.034152
At iterate 27 f = 2486.7 |proj g|= 0.043297
At iterate 28 f = 2486.7 |proj g|= 0.014902
At iterate 29 f = 2486.7 |proj g|= 0.016268
At iterate 30 f = 2486.7 |proj g|= 0.013669
At iterate 31 f = 2486.7 |proj g|= 0.021967
iterations 31
function evaluations 34
segments explored during Cauchy searches 54
BFGS updates skipped 0
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.0219674
final function value 2486.74
F = 2486.74
final value 2486.744525
converged
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~alpha_io + a_wl_io + bio_hum_cn + b_wl_io + dcatch_dlai_io +
dqcrit_io + dz0v_dh_io + f0_io + fd_io + g_area_io + g_root_io +
g_wood_io + gs_nvg_io + hw_sw_io + kaps_roth + knl_io + lai_max_io +
lai_min_io + lma_io + n_inorg_turnover + nmass_io + nr_io +
retran_l_io + retran_r_io + r_grow_io + rootd_ft_io + sigl_io +
sorp + tleaf_of_io + tlow_io + tupp_io + l_vg_soil
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -2315.205
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 2315.2 |proj g|= 1.9221
At iterate 1 f = 2314.2 |proj g|= 1.7645
At iterate 2 f = 2314.1 |proj g|= 0.74274
At iterate 3 f = 2313.9 |proj g|= 1.3088
At iterate 4 f = 2313.7 |proj g|= 1.9298
At iterate 5 f = 2313.5 |proj g|= 1.9378
At iterate 6 f = 2313.1 |proj g|= 1.9199
At iterate 7 f = 2312 |proj g|= 1.9412
At iterate 8 f = 2311.8 |proj g|= 1.9461
At iterate 9 f = 2310.4 |proj g|= 1.3481
At iterate 10 f = 2309.3 |proj g|= 1.2924
At iterate 11 f = 2308.7 |proj g|= 1.937
At iterate 12 f = 2308.3 |proj g|= 1.5008
At iterate 13 f = 2308.3 |proj g|= 1.216
At iterate 14 f = 2308.2 |proj g|= 1.8896
At iterate 15 f = 2308.1 |proj g|= 1.8935
At iterate 16 f = 2307.5 |proj g|= 1.9032
At iterate 17 f = 2307 |proj g|= 1.8988
At iterate 18 f = 2304.7 |proj g|= 1.0033
At iterate 19 f = 2303.6 |proj g|= 1.5714
At iterate 20 f = 2303.5 |proj g|= 1.917
At iterate 21 f = 2303.3 |proj g|= 1.9143
At iterate 22 f = 2303.2 |proj g|= 1.9086
At iterate 23 f = 2303.1 |proj g|= 1.8843
At iterate 24 f = 2302.9 |proj g|= 0.70255
At iterate 25 f = 2302.9 |proj g|= 0.68074
At iterate 26 f = 2302.9 |proj g|= 1.9109
At iterate 27 f = 2302.8 |proj g|= 1.7663
At iterate 28 f = 2302.8 |proj g|= 0.67488
At iterate 29 f = 2302.6 |proj g|= 1.0907
At iterate 30 f = 2302.3 |proj g|= 1.9095
At iterate 31 f = 2302.2 |proj g|= 1.9105
At iterate 32 f = 2302 |proj g|= 1.4736
At iterate 33 f = 2302 |proj g|= 1.4695
At iterate 34 f = 2301.8 |proj g|= 1.32
At iterate 35 f = 2301.7 |proj g|= 1.91
At iterate 36 f = 2301.7 |proj g|= 1.3149
At iterate 37 f = 2301.6 |proj g|= 0.83873
At iterate 38 f = 2301.5 |proj g|= 0.96389
At iterate 39 f = 2301.5 |proj g|= 1.3464
At iterate 40 f = 2301.4 |proj g|= 0.53138
At iterate 41 f = 2301.3 |proj g|= 1.8834
At iterate 42 f = 2301.3 |proj g|= 0.82272
At iterate 43 f = 2301.3 |proj g|= 0.75558
At iterate 44 f = 2301.3 |proj g|= 1.8803
At iterate 45 f = 2301.3 |proj g|= 1.8814
At iterate 46 f = 2301.3 |proj g|= 1.8812
At iterate 47 f = 2301.3 |proj g|= 1.881
At iterate 48 f = 2301.2 |proj g|= 1.3784
At iterate 49 f = 2301.2 |proj g|= 0.62431
At iterate 50 f = 2301.2 |proj g|= 1.9114
At iterate 51 f = 2301.2 |proj g|= 1.9123
At iterate 52 f = 2301.1 |proj g|= 1.9116
At iterate 53 f = 2301.1 |proj g|= 0.63241
At iterate 54 f = 2301.1 |proj g|= 0.28367
At iterate 55 f = 2301.1 |proj g|= 0.31082
At iterate 56 f = 2301.1 |proj g|= 0.62345
At iterate 57 f = 2301.1 |proj g|= 1.911
At iterate 58 f = 2301.1 |proj g|= 1.9119
At iterate 59 f = 2301.1 |proj g|= 1.9112
At iterate 60 f = 2301 |proj g|= 0.29321
At iterate 61 f = 2301 |proj g|= 0.31545
At iterate 62 f = 2301 |proj g|= 0.84934
At iterate 63 f = 2301 |proj g|= 1.8836
At iterate 64 f = 2301 |proj g|= 1.8821
At iterate 65 f = 2301 |proj g|= 0.58887
At iterate 66 f = 2301 |proj g|= 1.2003
At iterate 67 f = 2301 |proj g|= 0.43465
At iterate 68 f = 2301 |proj g|= 0.42034
At iterate 69 f = 2301 |proj g|= 0.60752
At iterate 70 f = 2301 |proj g|= 1.6659
At iterate 71 f = 2301 |proj g|= 1.8816
At iterate 72 f = 2300.9 |proj g|= 1.8816
At iterate 73 f = 2300.9 |proj g|= 1.8828
At iterate 74 f = 2300.9 |proj g|= 1.8811
At iterate 75 f = 2300.9 |proj g|= 0.40944
At iterate 76 f = 2300.9 |proj g|= 0.36834
At iterate 77 f = 2300.9 |proj g|= 0.74183
At iterate 78 f = 2300.9 |proj g|= 1.6322
At iterate 79 f = 2300.9 |proj g|= 1.914
At iterate 80 f = 2300.9 |proj g|= 1.9141
At iterate 81 f = 2300.9 |proj g|= 1.9155
At iterate 82 f = 2300.9 |proj g|= 1.1075
At iterate 83 f = 2300.8 |proj g|= 0.97412
At iterate 84 f = 2300.8 |proj g|= 0.96614
At iterate 85 f = 2300.8 |proj g|= 1.8816
At iterate 86 f = 2300.8 |proj g|= 1.8815
At iterate 87 f = 2300.8 |proj g|= 1.9143
At iterate 88 f = 2300.7 |proj g|= 0.61654
At iterate 89 f = 2300.7 |proj g|= 0.61993
At iterate 90 f = 2300.7 |proj g|= 0.61744
At iterate 91 f = 2300.7 |proj g|= 0.59652
At iterate 92 f = 2300.7 |proj g|= 0.47719
At iterate 93 f = 2300.7 |proj g|= 0.53494
At iterate 94 f = 2300.7 |proj g|= 0.53394
At iterate 95 f = 2300.7 |proj g|= 0.52042
At iterate 96 f = 2300.7 |proj g|= 0.60048
At iterate 97 f = 2300.7 |proj g|= 1.9142
At iterate 98 f = 2300.7 |proj g|= 1.9152
At iterate 99 f = 2300.6 |proj g|= 1.9161
At iterate 100 f = 2300.6 |proj g|= 1.9168
At iterate 101 f = 2300.6 |proj g|= 1.3656
final value 2300.595796
stopped after 101 iterations
Visualise the predicted outputs at the NROY points of the old design, and the new suggested design.
Interestingly, these aren’t perfectly within the original hard boundaries set by Andy, even though I’ve set those boundaries to be the +- 4 standard deviation threholds in the History Match. I suggest this is because there is model discrepancy, and that there is considerable wriggle room induced from emulator uncertainty.
In particular, it appears that vegetation carbon is difficult to keep high, and that many NROY proposed members have a fairly low vegetation carbon. This might need a discrepancy term, or adjusting in some other way. It certainly needs exploring, and a OAAT plot might give clues as to the parameters to choose.
How good are the four emulators that we’ve built? Are there biases? (there’s no real evidence of this)
It’s obviously hard to maintain a high vegetation carbon in particular. What parameter values might you choose to do this, and what might be the trade-offs you have to make?
First, what parameters affect vegetation carbon and how? How sure are we about that?
some values will be clipped
What might be the trade-offs for a high (or accurate) vegetation carbon? are they acceptable? Plot the oaat sensitivity of the other 3 outputs we’re calibrating on.
What do the emulators make of the design points which actually make Andy’s “hard boundary” criteria? If we leave them out, do they still place the output within the hard boundaries?
Do a leave-one-out cross validation of points inside the hard boundaries using the wave1 fits.
We see in the leave-one-out analysis that the emulator is consistently under-predicting the vegetation carbon (though the uncertainty estimate often covers the actual value).This suggests (1) that there isn’t really a huge problem with a model discrepancy (or at least that isn’t the only problem), and (2) the history matching is working as it should, and taking into account a not-great emulator.
Is a two-step emulator any better at emulating those crucial points which fall within aw’s hard boundaries? First, create a list of emulator fits.
It doesn’t appear that the two-step emulator (here plotted in red) is doing any better than the regular emulator.
Two other things I can think of to check: 1) how about using “multistart” to choose different starting conditions for optimising the emulators and 2) Using a flat prior for the mean function rather than a linear prior.
Loading required package: Matrix
Attaching package: ‘Matrix’
The following object is masked from ‘package:spam’:
det
Loaded glmnet 4.1-1
Input from Eddy Robertson:
“We can assume that the majority of vegetation carbon is stored in tree trunks so, the carbon density of trees is approximately (cVeg_lnd_mean / treeFrac_lnd_mean). Although I suppose it’s possible that in the PPE the shrubs have become tree-like, so it might be interesting to plot cVeg_lnd_mean versus (shrubFrac_lnd_mean + treeFrac_lnd_mean) as well. The first question is whether cVeg_lnd_mean increases linearly with treeFrac_lnd_mean. Then if this turned out to be interesting, we could produce a treeCVeg_lnd_mean output.”
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~1
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -324.7357
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 324.74 |proj g|= 1.8894
At iterate 1 f = 154.68 |proj g|= 2
At iterate 2 f = 144.37 |proj g|= 1.9477
At iterate 3 f = 140.12 |proj g|= 1.8272
At iterate 4 f = 137.65 |proj g|= 1.5752
At iterate 5 f = 135.35 |proj g|= 1.4243
At iterate 6 f = 134.81 |proj g|= 1.1667
At iterate 7 f = 134.52 |proj g|= 1.0068
At iterate 8 f = 134.47 |proj g|= 0.55475
At iterate 9 f = 134.44 |proj g|= 0.37207
At iterate 10 f = 134.43 |proj g|= 0.20987
At iterate 11 f = 134.42 |proj g|= 0.076796
At iterate 12 f = 134.42 |proj g|= 0.026753
At iterate 13 f = 134.42 |proj g|= 0.01759
At iterate 14 f = 134.42 |proj g|= 0.0060535
At iterate 15 f = 134.42 |proj g|= 0.0041787
At iterate 16 f = 134.42 |proj g|= 0.0021465
At iterate 17 f = 134.42 |proj g|= 0.00377
At iterate 18 f = 134.42 |proj g|= 0.001123
At iterate 19 f = 134.42 |proj g|= 0.00021622
iterations 19
function evaluations 20
segments explored during Cauchy searches 42
BFGS updates skipped 0
active bounds at final generalized Cauchy point 19
norm of the final projected gradient 0.000216221
final function value 134.421
F = 134.421
final value 134.420708
converged
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~1
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -1621.762
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 1621.8 |proj g|= 1.7706
At iterate 1 f = 1507.2 |proj g|= 1.4766
At iterate 2 f = 1423.3 |proj g|= 1.7768
At iterate 3 f = 1418.6 |proj g|= 1.4019
At iterate 4 f = 1415.9 |proj g|= 1.6446
At iterate 5 f = 1415.3 |proj g|= 1.7139
At iterate 6 f = 1415.1 |proj g|= 1.6554
At iterate 7 f = 1414.8 |proj g|= 1.4108
At iterate 8 f = 1414.7 |proj g|= 0.99795
At iterate 9 f = 1414.5 |proj g|= 0.73402
At iterate 10 f = 1414.5 |proj g|= 0.66567
At iterate 11 f = 1414.5 |proj g|= 0.13919
At iterate 12 f = 1414.5 |proj g|= 0.0858
At iterate 13 f = 1414.5 |proj g|= 0.034268
At iterate 14 f = 1414.5 |proj g|= 0.063708
At iterate 15 f = 1414.5 |proj g|= 0.010821
At iterate 16 f = 1414.5 |proj g|= 0.0075671
At iterate 17 f = 1414.5 |proj g|= 0.012372
At iterate 18 f = 1414.5 |proj g|= 0.0071278
At iterate 19 f = 1414.5 |proj g|= 0.0036723
iterations 19
function evaluations 21
segments explored during Cauchy searches 70
BFGS updates skipped 0
active bounds at final generalized Cauchy point 25
norm of the final projected gradient 0.00367225
final function value 1414.5
F = 1414.5
final value 1414.500017
converged
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~1
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -2800.645
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 2800.6 |proj g|= 1.8463
At iterate 1 f = 2616.7 |proj g|= 1.6218
At iterate 2 f = 2605.4 |proj g|= 1.553
At iterate 3 f = 2583.7 |proj g|= 1.6224
At iterate 4 f = 2577.2 |proj g|= 1.011
At iterate 5 f = 2573 |proj g|= 1.3781
At iterate 6 f = 2571.3 |proj g|= 1.308
At iterate 7 f = 2571.1 |proj g|= 0.9583
At iterate 8 f = 2571 |proj g|= 0.62272
At iterate 9 f = 2570.9 |proj g|= 0.44321
At iterate 10 f = 2570.9 |proj g|= 0.12021
At iterate 11 f = 2570.9 |proj g|= 0.055708
At iterate 12 f = 2570.9 |proj g|= 0.016351
At iterate 13 f = 2570.9 |proj g|= 0.0027779
iterations 13
function evaluations 16
segments explored during Cauchy searches 52
BFGS updates skipped 0
active bounds at final generalized Cauchy point 27
norm of the final projected gradient 0.00277786
final function value 2570.93
F = 2570.93
final value 2570.925355
converged
optimisation start
------------------
* estimation method : MLE
* optimisation method : BFGS
* analytical gradient : used
* trend model : ~1
* covariance model :
- type : matern5_2
- nugget : NO
- parameters lower bounds : 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
- parameters upper bounds : 2 2 1.999244 1.698141 1.994798 1.990622 1.998154 1.778477 2 2 2 1.998084 2 2 2 1.99598 2 1.989889 1.997502 1.989483 1.98997 1.987916 1.983347 1.988424 2 1.990947 1.995643 1.994232 1.992247 2 2 2
- best initial criterion value(s) : -2465.963
N = 32, M = 5 machine precision = 2.22045e-16
At X0, 0 variables are exactly at the bounds
At iterate 0 f= 2466 |proj g|= 1.7876
At iterate 1 f = 2357.6 |proj g|= 2
At iterate 2 f = 2342.2 |proj g|= 2
At iterate 3 f = 2334.4 |proj g|= 1.9202
At iterate 4 f = 2331.2 |proj g|= 1.8804
At iterate 5 f = 2325 |proj g|= 1.7272
At iterate 6 f = 2322 |proj g|= 1.6469
At iterate 7 f = 2320.5 |proj g|= 1.4532
At iterate 8 f = 2317.5 |proj g|= 1.4216
At iterate 9 f = 2317 |proj g|= 1.3362
At iterate 10 f = 2316.1 |proj g|= 1.2141
At iterate 11 f = 2315.5 |proj g|= 1.1842
At iterate 12 f = 2315.4 |proj g|= 0.93199
At iterate 13 f = 2315.4 |proj g|= 0.5059
At iterate 14 f = 2315.4 |proj g|= 0.38256
At iterate 15 f = 2315.4 |proj g|= 0.53325
At iterate 16 f = 2315.4 |proj g|= 0.18709
At iterate 17 f = 2315.4 |proj g|= 0.038636
At iterate 18 f = 2315.4 |proj g|= 0.02819
At iterate 19 f = 2315.4 |proj g|= 0.033483
At iterate 20 f = 2315.4 |proj g|= 0.015972
At iterate 21 f = 2315.4 |proj g|= 0.014807
At iterate 22 f = 2315.4 |proj g|= 0.0076813
At iterate 23 f = 2315.4 |proj g|= 0.0071707
iterations 23
function evaluations 27
segments explored during Cauchy searches 50
BFGS updates skipped 0
active bounds at final generalized Cauchy point 21
norm of the final projected gradient 0.00717074
final function value 2315.38
F = 2315.38
final value 2315.375074
converged
Is the range of carbon densities narrower in the “hard boundary” set?
For example, what effect has the AW “hard boundary” constraint produced on ranges of the other (marginal) outputs? Which constraint is doing the majority of the work?
A function to find the proportion of output space that is removed (or retained) when applying a constraint.
All constraints together take the normalized range of all the outputs to a little under a half of the original range, on average. Of course, four of the inputs are constrained here, which has quite a large impact.
We can see what impacts the constraints have individually
The NBP constraint has no marginal impact
Vegetation carbon (cVeg) has the strongest impact on everything else. This is no suprise, given that so few of the ensemble members seem to simulate carbon vegetation high enough.